Advancing the scale of synthetic biology via cross-species transfer of cellular functions enabled by iModulon engraftment.

Machine learning applied to large compendia of transcriptomic data has enabled the decomposition of bacterial transcriptomes to identify independently modulated sets of genes, such iModulons represent specific cellular functions. The identification of iModulons enables accurate identification of genes necessary and sufficient for cross-species transfer of cellular functions. We demonstrate cross-species transfer of: 1) the biotransformation of vanillate to protocatechuate, 2) a malonate catabolic pathway, 3) a catabolic pathway for 2,3-butanediol, and 4) an antimicrobial resistance to ampicillin found in multiple Pseudomonas species to Escherichia coli. iModulon-based engineering is a transformative strategy as it includes all genes comprising the transferred cellular function, including genes without functional annotation. Adaptive laboratory evolution was deployed to optimize the cellular function transferred, revealing mutations in the host. Combining big data analytics and laboratory evolution thus enhances the level of understanding of systems biology, and synthetic biology for strain design and development.

TCA cycle tailoring facilitates optimal growth of proton-pumping NADH dehydrogenase-dependent Escherichia coli.

IMPORTANCE: Energy generation pathways are a potential avenue for the development of novel antibiotics. However, bacteria possess remarkable resilience due to the compensatory pathways, which presents a challenge in this direction. NADH, the primary reducing equivalent, can transfer electrons to two distinct types of NADH dehydrogenases. Type I NADH dehydrogenase is an enzyme complex comprising multiple subunits and can generate proton motive force (PMF). Type II NADH dehydrogenase does not pump protons but plays a crucial role in maintaining the turnover of NAD+. To study the adaptive rewiring of energy metabolism, we evolved an Escherichia coli mutant lacking type II NADH dehydrogenase. We discovered that by modifying the flux through the tricarboxylic acid (TCA) cycle, E. coli could mitigate the growth impairment observed in the absence of type II NADH dehydrogenase. This research provides valuable insights into the intricate mechanisms employed by bacteria to compensate for disruptions in energy metabolism.

Laboratory evolution, transcriptomics, and modeling reveal mechanisms of paraquat tolerance.

Relationships between the genome, transcriptome, and metabolome underlie all evolved phenotypes. However, it has proved difficult to elucidate these relationships because of the high number of variables measured. A recently developed data analytic method for characterizing the transcriptome can simplify interpretation by grouping genes into independently modulated sets (iModulons). Here, we demonstrate how iModulons reveal deep understanding of the effects of causal mutations and metabolic rewiring. We use adaptive laboratory evolution to generate E. coli strains that tolerate high levels of the redox cycling compound paraquat, which produces reactive oxygen species (ROS). We combine resequencing, iModulons, and metabolic models to elucidate six interacting stress-tolerance mechanisms: (1) modification of transport, (2) activation of ROS stress responses, (3) use of ROS-sensitive iron regulation, (4) motility, (5) broad transcriptional reallocation toward growth, and (6) metabolic rewiring to decrease NADH production. This work thus demonstrates the power of iModulon knowledge mapping for evolution analysis.

Have you tried turning it off and on again? Oscillating selection to enhance fitness-landscape traversal in adaptive laboratory evolution experiments.

Adaptive Laboratory Evolution (ALE) is a powerful tool for engineering and understanding microbial physiology. ALE relies on the selection and enrichment of mutations that enable survival or faster growth under a selective condition imposed by the experimental setup. Phenotypic fitness landscapes are often underpinned by complex genotypes involving multiple genes, with combinatorial positive and negative effects on fitness. Such genotype relationships result in mutational fitness landscapes with multiple local fitness maxima and valleys. Traversing local maxima to find a global maximum often requires an individual or sub-population of cells to traverse fitness valleys. Traversing involves gaining mutations that are not adaptive for a given local maximum but are necessary to 'peak shift' to another local maximum, or eventually a global maximum. Despite these relatively well understood evolutionary principles, and the combinatorial genotypes that underlie most metabolic phenotypes, the majority of applied ALE experiments are conducted using constant selection pressures. The use of constant pressure can result in populations becoming trapped within local maxima, and often precludes the attainment of optimum phenotypes associated with global maxima. Here, we argue that oscillating selection pressures is an easily accessible mechanism for traversing fitness landscapes in ALE experiments, and provide theoretical and practical frameworks for implementation.

Adaptive evolution of a minimal organism with a synthetic genome.

The bacterial strain JCVI-syn3.0 stands as the first example of a living organism with a minimized synthetic genome, derived from the Mycoplasma mycoides genome and chemically synthesized in vitro. Here, we report the experimental evolution of a syn3.0- derived strain. Ten independent replicates were evolved for several hundred generations, leading to growth rate improvements of > 15%. Endpoint strains possessed an average of 8 mutations composed of indels and SNPs, with a pronounced C/G- > A/T transversion bias. Multiple genes were repeated mutational targets across the independent lineages, including phase variable lipoprotein activation, 5 distinct; nonsynonymous substitutions in the same membrane transporter protein, and inactivation of an uncharacterized gene. Transcriptomic analysis revealed an overall tradeoff reflected in upregulated ribosomal proteins and downregulated DNA and RNA related proteins during adaptation. This work establishes the suitability of synthetic, minimal strains for laboratory evolution, providing a means to optimize strain growth characteristics and elucidate gene functionality.

Adaptive laboratory evolution of Bacillus subtilis to overcome toxicity of lignocellulosic hydrolysate derived from Distiller's dried grains with solubles (DDGS).

Microbial tolerance to toxic compounds formed during biomass pretreatment is a significant challenge to produce bio-based products from lignocellulose cost effectively. Rational engineering can be problematic due to insufficient prerequisite knowledge of tolerance mechanisms. Therefore, adaptive laboratory evolution was applied to obtain 20 tolerant lineages of Bacillus subtilis strains able to utilize Distiller's Dried Grains with Solubles-derived (DDGS) hydrolysate. Evolved strains showed both improved growth performance and retained heterologous enzyme production using 100% hydrolysate-based medium, whereas growth of the starting strains was essentially absent. Whole-genome resequencing revealed that evolved isolates acquired mutations in the global regulator codY in 15 of the 19 sequenced isolates. Furthermore, mutations in genes related to oxidative stress (katA, perR) and flagella function appeared in both tolerance and control evolution experiments without toxic compounds. Overall, tolerance adaptive laboratory evolution yielded strains able to utilize DDGS-hydrolysate to produce enzymes and hence proved to be a valuable tool for the valorization of lignocellulose.

Laboratory evolution reveals general and specific tolerance mechanisms for commodity chemicals.

Although strain tolerance to high product concentrations is a barrier to the economically viable biomanufacturing of industrial chemicals, chemical tolerance mechanisms are often unknown. To reveal tolerance mechanisms, an automated platform was utilized to evolve Escherichia coli to grow optimally in the presence of 11 industrial chemicals (1,2-propanediol, 2,3-butanediol, glutarate, adipate, putrescine, hexamethylenediamine, butanol, isobutyrate, coumarate, octanoate, hexanoate), reaching tolerance at concentrations 60%-400% higher than initial toxic levels. Sequencing genomes of 223 isolates from 89 populations, reverse engineering, and cross-compound tolerance profiling were employed to uncover tolerance mechanisms. We show that: 1) cells are tolerized via frequent mutation of membrane transporters or cell wall-associated proteins (e.g., ProV, KgtP, SapB, NagA, NagC, MreB), transcription and translation machineries (e.g., RpoA, RpoB, RpoC, RpsA, RpsG, NusA, Rho), stress signaling proteins (e.g., RelA, SspA, SpoT, YobF), and for certain chemicals, regulators and enzymes in metabolism (e.g., MetJ, NadR, GudD, PurT); 2) osmotic stress plays a significant role in tolerance when chemical concentrations exceed a general threshold and mutated genes frequently overlap with those enabling chemical tolerance in membrane transporters and cell wall-associated proteins; 3) tolerization to a specific chemical generally improves tolerance to structurally similar compounds whereas a tradeoff can occur on dissimilar chemicals, and 4) using pre-tolerized starting isolates can hugely enhance the subsequent production of chemicals when a production pathway is inserted in many, but not all, evolved tolerized host strains, underpinning the need for evolving multiple parallel populations. Taken as a whole, this study provides a comprehensive genotype-phenotype map based on identified mutations and growth phenotypes for 223 chemical tolerant isolates.

Revealing oxidative pentose metabolism in new Pseudomonas putida isolates.

The Pseudomonas putida group in the Gammaproteobacteria has been intensively studied for bioremediation and plant growth promotion. Members of this group have recently emerged as promising hosts to convert intermediates derived from plant biomass to biofuels and biochemicals. However, most strains of P. putida cannot metabolize pentose sugars derived from hemicellulose. Here, we describe three isolates that provide a broader view of the pentose sugar catabolism in the P. putida group. One of these isolates clusters with the well-characterized P. alloputida KT2440 (Strain BP6); the second isolate clustered with plant growth-promoting strain P. putida W619 (Strain M2), while the third isolate represents a new species in the group (Strain BP8). Each of these isolates possessed homologous genes for oxidative xylose catabolism (xylDXA) and a potential xylonate transporter. Strain M2 grew on arabinose and had genes for oxidative arabinose catabolism (araDXA). A CRISPR interference (CRISPRi) system was developed for strain M2 and identified conditionally essential genes for xylose growth. A glucose dehydrogenase was found to be responsible for initial oxidation of xylose and arabinose in strain M2. These isolates have illuminated inherent diversity in pentose catabolism in the P. putida group and may provide alternative hosts for biomass conversion.

Understanding Functional Redundancy and Promiscuity of Multidrug Transporters in E. coli under Lipophilic Cation Stress.

Multidrug transporters (MDTs) are major contributors to microbial drug resistance and are further utilized for improving host phenotypes in biotechnological applications. Therefore, the identification of these MDTs and the understanding of their mechanisms of action in vivo are of great importance. However, their promiscuity and functional redundancy represent a major challenge towards their identification. Here, a multistep tolerance adaptive laboratory evolution (TALE) approach was leveraged to achieve this goal. Specifically, a wild-type E. coli K-12-MG1655 and its cognate knockout individual mutants ΔemrE, ΔtolC, and ΔacrB were evolved separately under increasing concentrations of two lipophilic cations, tetraphenylphosphonium (TPP+), and methyltriphenylphosphonium (MTPP+). The evolved strains showed a significant increase in MIC values of both cations and an apparent cross-cation resistance. Sequencing of all evolved mutants highlighted diverse mutational mechanisms that affect the activity of nine MDTs including acrB, mdtK, mdfA, acrE, emrD, tolC, acrA, mdtL, and mdtP. Besides regulatory mutations, several structural mutations were recognized in the proximal binding domain of acrB and the permeation pathways of both mdtK and mdfA. These details can aid in the rational design of MDT inhibitors to efficiently combat efflux-based drug resistance. Additionally, the TALE approach can be scaled to different microbes and molecules of medical and biotechnological relevance.

Experimental Evolution Reveals Unifying Systems-Level Adaptations but Diversity in Driving Genotypes.

Genotype-fitness maps of evolution have been well characterized for biological components, such as RNA and proteins, but remain less clear for systems-level properties, such as those of metabolic and transcriptional regulatory networks. Here, we take multi-omics measurements of 6 different E. coli strains throughout adaptive laboratory evolution (ALE) to maximal growth fitness. The results show the following: (i) convergence in most overall phenotypic measures across all strains, with the notable exception of divergence in NADPH production mechanisms; (ii) conserved transcriptomic adaptations, describing increased expression of growth promoting genes but decreased expression of stress response and structural components; (iii) four groups of regulatory trade-offs underlying the adjustment of transcriptome composition; and (iv) correlates that link causal mutations to systems-level adaptations, including mutation-pathway flux correlates and mutation-transcriptome composition correlates. We thus show that fitness landscapes for ALE can be described with two layers of causation: one based on system-level properties (continuous variables) and the other based on mutations (discrete variables). IMPORTANCE Understanding the mechanisms of microbial adaptation will help combat the evolution of drug-resistant microbes and enable predictive genome design. Although experimental evolution allows us to identify the causal mutations underlying microbial adaptation, it remains unclear how causal mutations enable increased fitness and is often explained in terms of individual components (i.e., enzyme rate) as opposed to biological systems (i.e., pathways). Here, we find that causal mutations in E. coli are linked to systems-level changes in NADPH balance and expression of stress response genes. These systems-level adaptation patterns are conserved across diverse E. coli strains and thus identify cofactor balance and proteome reallocation as dominant constraints governing microbial adaptation.

Laboratory evolution of synthetic electron transport system variants reveals a larger metabolic respiratory system and its plasticity.

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. There is a detailed understanding of the structural and biochemical features of respiratory enzymes; however, a holistic examination of the system and its plasticity is lacking. Here we generate four strains of Escherichia coli harboring unbranched ETS that pump 1, 2, 3, or 4 proton(s) per electron and characterized them using a combination of synergistic methods (adaptive laboratory evolution, multi-omic analyses, and computation of proteome allocation). We report that: (a) all four ETS variants evolve to a similar optimized growth rate, and (b) the laboratory evolutions generate specific rewiring of major energy-generating pathways, coupled to the ETS, to optimize ATP production capability. We thus define an Aero-Type System (ATS), which is a generalization of the aerobic bioenergetics and is a metabolic systems biology description of respiration and its inherent plasticity.

Engineering Pseudomonas putida for improved utilization of syringyl aromatics.

Lignin is a largely untapped source for the bioproduction of value-added chemicals. Pseudomonas putida KT2440 has emerged as a strong candidate for bioprocessing of lignin feedstocks due to its resistance to several industrial solvents, broad metabolic capabilities, and genetic amenability. Here we demonstrate the engineering of P. putida for the ability to metabolize syringic acid, one of the major products that comes from the breakdown of the syringyl component of lignin. The rational design was first applied for the construction of strain Sy-1 by overexpressing a native vanillate demethylase. Subsequent adaptive laboratory evolution (ALE) led to the generation of mutations that achieved robust growth on syringic acid as a sole carbon source. The best mutant showed a 30% increase in the growth rate over the original engineered strain. Genomic sequencing revealed multiple mutations repeated in separate evolved replicates. Reverse engineering of mutations identified in agmR, gbdR, fleQ, and the intergenic region of gstB and yadG into the parental strain recaptured the improved growth of the evolved strains to varied extent. These findings thus reveal the ability of P. putida to utilize lignin more fully as a feedstock and make it a more economically viable chassis for chemical production.

Membrane transporter identification and modulation via adaptive laboratory evolution.

Membrane transport proteins are potential targets for medical and biotechnological applications. However, more than 30% of reported membrane transporter families are either poorly characterized or lack adequate functional annotation. Here, adaptive laboratory evolution was leveraged to identify membrane transporters for a set of four amino acids as well as specific mutations that modulate the activities of these transporters. Specifically, Escherichia coli was adaptively evolved under increasing concentrations of L-histidine, L-phenylalanine, L-threonine, and L-methionine separately with multiple replicate evolutions. Evolved populations and isolated clones displayed growth rates comparable to the unstressed ancestral strain at elevated concentrations (four-to six-fold increases) of the targeted amino acids. Whole genome sequencing of the evolved strains revealed a diverse number of key mutations, including SNPs, small deletions, and copy number variants targeting the transporters leuE for histidine, yddG for phenylalanine, yedA for methionine, and brnQ and rhtC for threonine. Reverse engineering of the mutations in the ancestral strain established mutation causality of the specific mutations for the tolerant phenotypes. The functional roles of yedA and brnQ in the transport of methionine and threonine, respectively, are novel assignments and their functional roles were validated using a flow cytometry cellular accumulation assay. To demonstrate how the identified transporters can be leveraged for production, an L-phenylalanine overproduction strain was shown to be a superior producer when the identified yddG exporter was overexpressed. Overall, the results revealed the striking efficiency of laboratory evolution to identify transporters and specific mutational mechanisms to modulate their activities, thereby demonstrating promising applicability in transporter discovery efforts and strain engineering.

Machine-learning from Pseudomonas putida KT2440 transcriptomes reveals its transcriptional regulatory network.

Bacterial gene expression is orchestrated by numerous transcription factors (TFs). Elucidating how gene expression is regulated is fundamental to understanding bacterial physiology and engineering it for practical use. In this study, a machine-learning approach was applied to uncover the genome-scale transcriptional regulatory network (TRN) in Pseudomonas putida KT2440, an important organism for bioproduction. We performed independent component analysis of a compendium of 321 high-quality gene expression profiles, which were previously published or newly generated in this study. We identified 84 groups of independently modulated genes (iModulons) that explain 75.7% of the total variance in the compendium. With these iModulons, we (i) expand our understanding of the regulatory functions of 39 iModulon associated TFs (e.g., HexR, Zur) by systematic comparison with 1993 previously reported TF-gene interactions; (ii) outline transcriptional changes after the transition from the exponential growth to stationary phases; (iii) capture group of genes required for utilizing diverse carbon sources and increased stationary response with slower growth rates; (iv) unveil multiple evolutionary strategies of transcriptome reallocation to achieve fast growth rates; and (v) define an osmotic stimulon, which includes the Type VI secretion system, as coordination of multiple iModulon activity changes. Taken together, this study provides the first quantitative genome-scale TRN for P. putida KT2440 and a basis for a comprehensive understanding of its complex transcriptome changes in a variety of physiological states.

Selection for Cell Yield Does Not Reduce Overflow Metabolism in Escherichia coli.

Overflow metabolism is ubiquitous in nature, and it is often considered inefficient because it leads to a relatively low biomass yield per consumed carbon. This metabolic strategy has been described as advantageous because it supports high growth rates during nutrient competition. Here, we experimentally evolved bacteria without nutrient competition by repeatedly growing and mixing millions of parallel batch cultures of Escherichia coli. Each culture originated from a water-in-oil emulsion droplet seeded with a single cell. Unexpectedly we found that overflow metabolism (acetate production) did not change. Instead, the numerical cell yield during the consumption of the accumulated acetate increased as a consequence of a reduction in cell size. Our experiments and a mathematical model show that fast growth and overflow metabolism, followed by the consumption of the overflow metabolite, can lead to a higher numerical cell yield and therefore a higher fitness compared with full respiration of the substrate. This provides an evolutionary scenario where overflow metabolism can be favorable even in the absence of nutrient competition.

Escherichia coli Data-Driven Strain Design Using Aggregated Adaptive Laboratory Evolution Mutational Data.

Microbes are being engineered for an increasingly large and diverse set of applications. However, the designing of microbial genomes remains challenging due to the general complexity of biological systems. Adaptive Laboratory Evolution (ALE) leverages nature's problem-solving processes to generate optimized genotypes currently inaccessible to rational methods. The large amount of public ALE data now represents a new opportunity for data-driven strain design. This study describes how novel strain designs, or genome sequences not yet observed in ALE experiments or published designs, can be extracted from aggregated ALE data and demonstrates this by designing, building, and testing three novel Escherichia coli strains with fitnesses comparable to ALE mutants. These designs were achieved through a meta-analysis of aggregated ALE mutations data (63 Escherichia coli K-12 MG1655 based ALE experiments, described by 93 unique environmental conditions, 357 independent evolutions, and 13 957 observed mutations), which additionally revealed global ALE mutation trends that inform on ALE-derived strain design principles. Such informative trends anticipate ALE-derived strain designs as largely gene-centric, as opposed to noncoding, and composed of a relatively small number of beneficial variants (approximately 6). These results demonstrate how strain design efforts can be enhanced by the meta-analysis of aggregated ALE data.

Machine Learning of Bacterial Transcriptomes Reveals Responses Underlying Differential Antibiotic Susceptibility.

In vitro antibiotic susceptibility testing often fails to accurately predict in vivo drug efficacies, in part due to differences in the molecular composition between standardized bacteriologic media and physiological environments within the body. Here, we investigate the interrelationship between antibiotic susceptibility and medium composition in Escherichia coli K-12 MG1655 as contextualized through machine learning of transcriptomics data. Application of independent component analysis, a signal separation algorithm, shows that complex phenotypic changes induced by environmental conditions or antibiotic treatment are directly traced to the action of a few key transcriptional regulators, including RpoS, Fur, and Fnr. Integrating machine learning results with biochemical knowledge of transcription factor activation reveals medium-dependent shifts in respiration and iron availability that drive differential antibiotic susceptibility. By extension, the data generation and data analytics workflow used here can interrogate the regulatory state of a pathogen under any measured condition and can be applied to any strain or organism for which sufficient transcriptomics data are available. IMPORTANCE Antibiotic resistance is an imminent threat to global health. Patient treatment regimens are often selected based on results from standardized antibiotic susceptibility testing (AST) in the clinical microbiology lab, but these in vitro tests frequently misclassify drug effectiveness due to their poor resemblance to actual host conditions. Prior attempts to understand the combined effects of drugs and media on antibiotic efficacy have focused on physiological measurements but have not linked treatment outcomes to transcriptional responses on a systems level. Here, application of machine learning to transcriptomics data identified medium-dependent responses in key regulators of bacterial iron uptake and respiratory activity. The analytical workflow presented here is scalable to additional organisms and conditions and could be used to improve clinical AST by identifying the key regulatory factors dictating antibiotic susceptibility.

Genome-scale metabolic modeling reveals key features of a minimal gene set.

Mesoplasma florum, a fast-growing near-minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein-coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species-specific constraints to produce the functional iJL208 genome-scale model (GEM) of metabolism. Genome-wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI-syn3.0 and its parent JCVI-syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping-stone toward model-driven whole-genome engineering.

Environmental conditions dictate differential evolution of vancomycin resistance in Staphylococcus aureus.

While microbiological resistance to vancomycin in Staphylococcus aureus is rare, clinical vancomycin treatment failures are common, and methicillin-resistant S. aureus (MRSA) strains isolated from patients after prolonged vancomycin treatment failure remain susceptible. Adaptive laboratory evolution was utilized to uncover mutational mechanisms associated with MRSA vancomycin resistance in a physiological medium as well as a bacteriological medium used in clinical susceptibility testing. Sequencing of resistant clones revealed shared and media-specific mutational outcomes, with an overlap in cell wall regulons (walKRyycHI, vraSRT). Evolved strains displayed similar properties to resistant clinical isolates in their genetic and phenotypic traits. Importantly, resistant phenotypes that developed in physiological media did not translate into resistance in bacteriological media. Further, a bacteriological media-specific mechanism for vancomycin resistance associated with a mutated mprF was confirmed. This study bridges the gap between the understanding of clinical and microbiological vancomycin resistance in S. aureus and expands the number of allelic variants (18 ± 4 mutations for the top 5 mutated genes) that result in vancomycin resistance phenotypes.

Compensatory evolution of Pseudomonas aeruginosa's slow growth phenotype suggests mechanisms of adaptation in cystic fibrosis.

Long-term infection of the airways of cystic fibrosis patients with Pseudomonas aeruginosa is often accompanied by a reduction in bacterial growth rate. This reduction has been hypothesised to increase within-patient fitness and overall persistence of the pathogen. Here, we apply adaptive laboratory evolution to revert the slow growth phenotype of P. aeruginosa clinical strains back to a high growth rate. We identify several evolutionary trajectories and mechanisms leading to fast growth caused by transcriptional and mutational changes, which depend on the stage of adaptation of the strain. Return to high growth rate increases antibiotic susceptibility, which is only partially dependent on reversion of mutations or changes in the transcriptional profile of genes known to be linked to antibiotic resistance. We propose that similar mechanisms and evolutionary trajectories, in reverse direction, may be involved in pathogen adaptation and the establishment of chronic infections in the antibiotic-treated airways of cystic fibrosis patients.